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ENERGY DEPENDENCE OF ELECTRON-INDUCED 
RADIATION DAMAGE IN TUNGSTEN 
by James A. DiCarlo and James T. Stanley* 

Lewis Research Center 

SUMMARY 

The energy dependence of low dose damage production in commercial and high purity 
polycrystalline tungsten wires was studied near 350 K with 1. 6 to 2.4 MeV electrons. 
From resistivity measurements at 291 K the threshold energy for the onset of observable 
damage was determined as 50±2 eV. An "effective" threshold of 52±2 eV was also 
determined by directly fitting the energy dependence of the damage rates to theo- 
retical displacement cross sections calculated from step-function displacement prob- 
abilities. A decrease of two orders of magnitude in impurity content reduced damage 
rates by about a factor of two but did not affect threshold. These results combined 
with current defect recovery models for tungsten, low temperature threshold data, 
and computer-calculated bcc damage theory suggest: (1) Observed damage consisted of 
equal concentrations of vacancies and impurity -trapped Stage I free interstitials. 

(2) Across Stage II (100 to 600 K) onset threshold should be within 50±2 eV. (3) Minimum 
recoil energy required for free interstitial production near 0 K is 53±5 eV. (4) Thresh- 
old has little dependence on crystal direction. 

An empirical method is presented for predicting threshold energies in the bcc tran- 
sition metals by assuming the directional dependence of threshold is directly proportional 
to that of Young's modulus. By the use of one universal proportionality constant 
(1.2 x10” d (eV)(cm^)/N), thresholds for a number of metals and directions are calculated 
and shown to have significantly better agreement with experiment than the best available 
theoretical estimates. 


* Associate Professor of Engineering, Arizona State University, Tempe, Arizona. 


INTRODUCTION 


The present study was undertaken to determine the effects of a Stage II irradiation 
temperature on the threshold energy for displacement and the production of electron- 
induced damage in tungsten. Stage II for tungsten is that defect recovery stage from 100 
to 600 K which has been attributed to the release of impurity -trapped interstitials that 
first become mobile during Stage I (refs. 1 and 2). Although only a fraction of the total 
damage is retained, irradiations within Stage II can provide useful information concern- 
ing the effects of temperature and purity on the threshold energy and the damage rate. 
From these results further information might then be deduced concerning the recoil en- 
ergy required for free interstitial production and the dependence of threshold energy on 
recoil direction. The Stage II irradiations of the present work were performed at ~350 K 
not only for experimental convenience but also because the knowledge obtained near room 
temperature is more technically useful than that obtained at lower temperatures. 


THEORY 

The threshold energy for displacement is the minimum recoil energy which an 
atom in a perfect lattice must receive in order to produce an atomic displacement. The 
simplest manner in which it can be determined is by measuring the electron energy e ^ 
at which the first observable irradiation -induced damage is produced. Then, from en- 
ergy transfer considerations, the threshold energy for the onset of damage is given 
by 


Ej = E 


max 


(* 2 > 
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2 

E max^ _ 2(e + 2mc ) e ^ 

r 9 

Mc z 

Here E is the maximum recoil energy that the struck atom can receive from an 
r 2 2 

electron of energy e, and me and Me are the rest mass energy of the electron and 
atom, respectively. 

At irradiation temperatures within the lower portion of Stage I, there are situations 
when the damage onset method cannot be employed for an accurate threshold determina- 
tion. This is due primarily to the existence of subthreshold displacements caused by 
extrinsic sources (ref. 3). In these situations the usual approach is to assume that the 
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probability for displacement P(E r ) is a step function, that is, P = 0 for E r < E*j and 
P = 1 for E r > E^. The ’’effective" threshold energy E*j can then be determined by 
measuring the energy dependence of the damage rate and finding the ’’best" theoretical 
displacement cross section which obeys the relation 


2|(0=P F a d (6, E*) 


( 3 ) 


where dp/d$, the damage rate, is the change in specimen resistivity p per unit change 
in electron dose «i>; p F is the specific resistivity per interstitial -vacancy (Frenkel) 
pair; and cr d {e,E d ) is the displacement cross section calculated from a step function 
P(E r ). Here it is assumed that all displaced atoms are retained and observed and that 
p F is independent of interstitial -vacancy separation. For many metals in which both 
thresholds have been determined, E*j has been found to be equal to or slightly larger than 
E^ by only a few electron volts (ref. 3). 

Because equation (3) assumes the retention of all damage, it cannot describe the en- 
ergy and dose dependence of the damage rate for irradiation temperatures above Stage I 
defect recovery. The proper Stage n relation which has had experimental verification 
(refs. 4 and 5) was originally given by Walker (ref. 6); namely, 
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Here f(e) is the fraction of displaced atoms or interstitials created far enough from their 
own vacancies to be uncorrelated; p is the resistivity per Frenkel pair when the inter- 
stitial member is trapped at an impurity atom; and N u is the total concentration of im- 
purity atoms capable of capturing interstitials at the irradiation temperature. These im- 
purity atoms are considered to be "unsaturable"; that is, their trapping efficiency is un- 
altered after capturing one or more interstitials. 

The right side of equation (4) can be understood as the product of three probabilities: 
the probability per unit dose of producing a displaced atom - cr d (e,E d ); the probability 
that this displacement is an uncorrelated interstitial - f(e); and the probability that this 
interstitial will be captured by an unsaturable trap before encountering a vacancy - 
Jl + (2a d $/N u )J "l/^. The physical basis for the latter two probabilities centers on the 
model that at a Stage II temperature all irradiation-induced interstitials have sufficient 
thermal energy to migrate immediately after creation. Most perform correlated motion 
toward their own vacancies where they disappear from experimental detection. However, 
a certain fraction f of the interstitials are uncorrelated or "free" interstitial atoms 
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(FIA) which, because they wander randomly in the lattice, can be trapped by an impurity 
before encountering a vacancy. Thus the impurities retain the FIA and their vacancies 
for observation. Walker assumes that the probability for FLA to be trapped is given by 
N u /(N u + v) where v is the total vacancy concentration at the time the FIA are wandering 
around the lattice. Initially N u » v so that essentially no FIA reach vacancies and 
resistivity increases linearly with dose. As v increases, however, more and more FIA 
reach and annihilate vacancies before becoming trapped, resulting in an observable re- 
duction in dp/d4>. Eventually v » N u and dp/d$ = 0. Since v = a^t>, competition be- 
tween v and N u for FIA capture is contained mathematically in the denominator of 
equation (4). 

Determination of with equation (4) is complicated by the parameters N u and f. 
To minimize their effects, the following simplifications were employed in the present 
study. First, because N u could not be accurately estimated, damage rates were meas- 
ured only in the low dose region where N u » v and dp/d4> is independent of 4>. Under 
these conditions equation (4) reduces to 


PW =f(£)p*<J d (£,E*) (5) 

Ld4> J n 

Secondly, because f(e) should increase with energy in an unknown manner, the approach 
taken here is to define E*j for Stage II as that effective threshold energy which satisfies 
equation (5) under the assumption that f(e) = f = constant. That is-, for a Stage II irra- 
diation temperature T, the damage rate is considered to be directly proportional to the 

j|£ 

displacement cross section so that E,(T) can be determined by the same methods dis- 

d * 

cussed with equation (3). It should be pointed out that the E^ calculated in this manner 
may differ from the value calculated with the full theory (eq. (4)). However, as the ex- 
perimental results will show, this approach does yield an E^fT) essentially equal to 
E^j(T) and thus, at least for tungsten, is another and perhaps more convenient method for 
determining the Stage II threshold energy for displacement. 


EXPERIMENTAL PROCEDURE 

The specimens primarily employed for the threshold determinations were commer- 
cial purity 2 -mil -diameter polycrystalline wires supplied by General Electric. Resis- 
tivity ratio measurements yielded values of about 10, indicating a purity of ~99. 5 percent 
for this material (ref. 7). A high purity 5-mil-diameter polycrystalline wire was used to 
check both the commercial purity results and previous low temperature data; that is, the 
high purity wire was identical to that employed by DiCar lo, Snead, and Goland (ref. 8) 
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for radiation damage studies near 20 K. The residual resistivity ratio of this wire was 
estimated to be greater than 1000, indicating a purity of ~99. 995 percent and thus an im- 
purity concentration some two orders of magnitude smaller than the commercial wire. 

Resistivity changes caused by the electron damage were measured at 291 K using a 
specimen and dummy in a standard bridge circuit to eliminate thermal effects. Each 
wire was held at one end by a copper clamp and at the other end by a mica-insulated 
clamp, resulting in an exposed length of about 1. 2 centimeters. Voltage probes con- 
sisted of 1-mil tungsten wires spot -welded 1 centimeter apart. Because the ratio of 
measurement length to cross-sectional area could not be made identical for the two 
wires, there existed the possibility of a slight thermal resistivity error if consecutive 
measurements were not made at exactly the same temperature. This type of error was 
minimized in two ways. First, the preirradiation currents through the specimen and 
dummy were separately adjusted until the voltage drops across each wire were identical. 
Theoretically it can be shown that if these preirradiation currents were maintained dur- 
ing subsequent measurements, errors due to nonrepetitive temperatures should not exist. 
Experimentally, however, this method was not completely adequate since resistivity 
deviations could be observed if the temperature variations were allowed to exceed 1 K. 
Therefore a temperature controller was introduced to repeat temperatures to within \ 

±0. 2 K as measured by a Chromel-Alumel thermocouple. The controller, a PAR model 
152, balanced the input current to heaters within the specimen holder against the cooling 
power of a thermoelectric cooler. Resistivity and current ratio measurements were > 
made using a 6 -dial Rubicon double potentiometer with a voltage resolution of 0. 01 micro- 
volt. Since the currents through the 2- and 5-mil specimens were maintained at 40 and 
100 milliamperes, respectively, the total measurement system had a resolution on the 
order of 10 - ^ ohm -centimeter in defect resistivity. 

Electron irradiations were performed in the vacuum of the Lewis Dynamitron accel- 
erator. Near beam uniformity was achieved by magnetically expanding the beam and al- 
lowing about 10 percent or less of its center portion to pass through an aperature imme- 
diately in front of the specimen. The consistency of this uniformity was continually 
checked by monitoring the specimen’s voltage rise under identical beam current condi- 
tions. It was determined both theoretically and experimentally that beam heating caused 
the specimens to rise to an average temperature of about 350 K. This temperature was 
produced by an operating beam current density of ~50 microamperes per square centi- 
meter measured by a Faraday cup directly behind the specimens. Systematic and ran- 
dom errors in the integrated dose measurements were estimated at less than 5 and 1 per- 
cent, respectively. 

Electron energy was measured by the divider -board current of the Lewis Dynamitron 
accelerator and the 1968 calibration results of Stevenson and Gauntner (ref. 9). A 
Be(y,n) calibration at 1. 662 MeV indicated that since this latter study the resistance of 
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the divider board had not changed more than 0. 1 percent (~1. 7 keV). From a considera- 
tion of this result and other possible sources of error, each energy measurement of the 
present study was estimated to have a total absolute error less than 10 keV. 

For an accurate determination of damage rate, the measured values of dp/d<3> had 
to be corrected for two energy -dependent deviations caused by the finite diameter of the 
wire specimens. First, due to electron energy losses within the wires, the displace- 
ment cross section decreases with depth into the material, giving rise to position- 
dependent resistivity changes. According to Sosin (ref. 10), this energy loss effect can 
be averaged by assuming that the resistivity was measured at an average electron energy 
e given by 


e = e- - — cud (6) 

1 3t r 

where Cj is the incident electron energy, a is the specific energy loss per unit dis- 
tance, and d is the diameter of each wire. In the present work the a's calculated by 
Berger and Seltzer (ref. 11) were employed. 

The second energy -dependent correction is required because the electrons undergo 
multiple scatterings within each specimen, resulting in an average electron path length 
greater than the specimen thickness parallel to the beam direction. Because an in- 
creased path length implies a larger effective electron flux dependent on energy and posi- 
tion, the measured resistivity damage rate (dp/d«f>) m is an energy- and thickness - 
dependent fraction of (dp/d^)^, the "true” or zero diameter damage rate. Again accord- 
ing to Sosin (ref. 10), 


m + (7 

\d#/ 0 \d$/ m \ 4 / 

The parameter a can be determined from the average path length increase A in a foil 
of thickness t; that is, 


a = ^i) 

t 2 

As discussed by Oen (ref. 12) there has been some question as to the correct theoretical 
expression for A, often referred to as the Yang correction. Following Oen’s suggestion 
the present work employs the following relation for a: 
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( 8 ) 


a = m = [ e w ( j 2 

T 2 4t 

where t (= 7rd/4) is the average thickness of each wire and is the empirically cor- 
rected mean scattering angle contained in Moliere multiple scattering theory (ref. 13). 


EXPERIMENTAL RESULTS 


The general measurement procedure for determining damage rate as a function of 
energy was as follows. Because preirradiation handling may have introduced defects 
which anneal at the irradiation temperatures (ref. 7), each specimen was initially beam- 
annealed for a few minutes at an average temperature of ~420 K. After the anneal, the 
specimen was irradiated at constant energy until linearity of the resistivity production 
curve was established and an accurate value of the slope could be determined. A new en- 
ergy was then chosen and the same specimen was reirradiated without any intervening 
anneal. After slope measurement the aforementioned process was repeated until the 

_ O 

total irradiation-induced resistivity reached an arbitrary value of ~10 ohm -centimeter 
(~2xi0 -5 defect concentration). Typical resistivity production curves are given in fig- 
ure 1 for one of the commercial specimens. Table I lists the slope and standard devia- 
tion of each production curve at the incident energy and the corresponding maximum 
recoil energy E™ 3 ^ given by equation (2). The deviations in slope were attributed pri- 
marily to slight variations in beam uniformity, except for the first commercial specimen 
where no temperature controller was employed. 

From table I it can be seen that 350 K radiation damage in tungsten does not meas- 
urably commence in the commercial wires until electrons attain an energy somewhere 
between 1. 58 and ~1. 66 MeV. This latter value was inferred from a worst case extrap- 
olation of the finite damage rate measured at 1. 70±0. 01 MeV. Therefore, using equa- 
tions (1) and (2) one can conclude that near 350 K the onset threshold energy for displace- 
ment in tungsten lies between 48 and 52 eV; that is, E^(350 K) = 50±2 eV. 

In order to compute the effective threshold energy, the energy loss and path length 
corrections given by equations (6) to (8) were applied to the data of table I. The results 
are plotted in figure 2. The commercial data points are indicated by open blocks and the 
limited high purity data by solid blocks. From equation (5) and the tabulated displace- 
ment cross section data of Burke, Grossbard, and Lowe (ref. 14), it was determined that 
the best fit u^(e,E^) for the commercial results gave E*(350 K) = 52±2 eV and 
f Q p* = (2. I±0.4)xl0 - ^ ohm -centimeter. The energy dependence of the product 
f 0 p*a d (e,52 eV) is shown by the solid line of figure 2. The best fit cross sections were 
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determined only with the data above 1. 9 MeV since deviations in threshold energy due to 
lattice structure, subthreshold displacements, and specimen thickness may be present in 
the near -threshold data (ref. 3). However, as can be seen, the cr^ based on a 52 -eV 
step function fits nearly all the low energy data. Thus, even for irradiation temper- 
atures above Stage I, the threshold determination method which assumes the (low dose) 
damage rate to be directly proportional to a theoretical a^(E^) can yield an E d (T) essen- 
tially equal to E^(T). 

As shown by the dashed line of figure 2, the best fit curve for the high purity damage 
^ + 4 

rates yields E^ = 54 eV and f Q p =1.3x10 ohm -centimeter. Although agreement with 
the commercial threshold is good, the fact that there are only two points makes any E*j 
assignment for this material somewhat tenuous. Also, the validity of applying the path 
length correction to the high purity data is questionable since for 5 -mil tungsten 
(~60°) is well above the maximum Moliere scattering angle (30°). However, assuming 
that only the magnitude and not the energy dependence of the path length correction was 
in error, it was determined that doubling or even eliminating the correction has negli- 
gible effect on E^ for either specimen type but only alters f p . It should also be 
pointed out that, no matter what its absolute value may be, the path length correction 
should have a greater reducing effect on the measured damage rates of the larger diam- 
eter high purity specimen than on the measured rates of the commercial specimens. 

Thus the smaller table I damage rates for the high purity specimen are real and cannot 
be accounted for by the path length correction. 


DISCUSSION 

The threshold energy and damage rate results of the present study have both basic 
and applied implications which can only be appreciated in the light of radiation damage 
theory and other experimental data. For instance, by comparing the 350 K threshold en- 
ergies with similar low temperature results, one may obtain fundamental knowledge con- 
cerning the effects of temperature, recoil energy, and recoil direction on the production 
and deposition of damage. Although this knowledge would primarily concern tungsten, 
inferences for other bcc metals can result. In addition, there is the practical aspect of 
the damage rate results, namely, the capability of estimating total irradiation-induced 
defect concentrations when tungsten of a certain purity is subjected to electron bombard- 
ment at Stage n temperatures. The discussion that follows will consider each of the 
aforementioned topics in more detail. 
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THRESHOLD IN TUNGSTEN 
Theoretical Effects of Temperature 

Before discussing the experimental temperature dependence of threshold in tungsten, 
one should first consider the possible theoretical effects of temperature on the displace- 
ment and retention of an atom in a real crystal. The many aspects of this subject have 
been adequately covered in a recent review article by Sosin and Bauer (ref. 3). Hence 
we will briefly re-examine only that theory pertinent to the present experiment, concen- 
trating on such topics as displacement theory with and without thermal vibration and the 
mechanisms currently invoked for low temperature defect recovery in tungsten. 

Disp lacement theory . - Consider defect production in a perfect lattice at 0 K. Be- 
cause the structure surrounding each atom in the lattice is discrete, the threshold energy 
for displacement is a function of the atom's initial recoil direction; that is, 

= E^hkZ) where the indices h, k, and l describe the recoil direction in terms of 
the principal crystal axes. Since E^(0 K) is the minimum recoil energy required to pro- 
duce a displacement, it follows that E^hkZ) > E^(0 K). Computer studies on both fee 
(ref. 15) and bcc (ref. 16) lattices have found that in contrast to the vacancy the intersti- 
tial is not created at the collision site but is usually the result of a replacement collision 
chain. That is, the initial recoil atom first displaces one of its neighbors by transfer- 
ring a large fraction of E r and then replaces this neighbor in its original position. This 
sequence of energy transfer and replacement continues from neighbor to neighbor until a 
collision occurs in which the energy transferred is insufficient to displace an atom from 
its original site. The energy loss per replacement was found to vary with crystal direc- 
tion of the chain, being smallest for the closest -packed direction. Thus the final out- 
come of the electron collision is an atom in an interstitial position separated from a va- 
cant lattice site by a distance which is dependent on the initial recoil energy, the initial 
recoil direction, and the crystal direction in which the replacement chain occurred. 

If E r , the maximum recoil energy available from the incident electrons (eq. (2)), 
equals E£j(0 K), the interstitial will be created as close to the vacancy as possible with- 
out recombining. Therefore E^(0 K) is the threshold energy required for the production 
of a Frenkel pair characterized by minimum separation between the vacancy and the 
"close-pair" interstitial atom (CPIA). On the average this close-pair separation should 
increase with E^ 13 ^ - E^j(0 K). Eventually at some E^ 1 ^ equal to E ( j(FIA) a few inter- 
stitials will attain a separation large enough to be outside the attractive influence of their 
own vacancies, thereby becoming uncorrelated or FIA. 

Introducing thermal vibration into a perfect lattice can affect the two basic aspects 
of damage production, that is, the displacement of the initial recoil atom and the reten- 
tion of the resulting interstitial. For example, temperature can alter the displacement 
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process by contributing thermal momentum to the collision momentum of the recoil 
atom, thereby lowering the onset threshold energy (ref. 3). If the collision and thermal 
momentum happen to be parallel at impact, then 

E d (T) “ E d (0 K) “ 2 ( E d E v) 1/2 (9) 

where E y = kT is the thermal energy of the struck atom. For = 50 eV and 
T = 350 K, E^(0 K) - E^(T) = 2. 5 eV. If the thermal momentum of the atoms surround- 
ing the struck atom happen to be 180° out of phase with that of the atom, another 2. 5 eV 
energy gain may result. Thus, due to this temperature -enhanced displacement effect 
(TEDE), E^ measured at 350 K in tungsten may be as much as 5 eV lower than that 
measured near 0 K. However, Et should be unaffected by TEDE since it is determined 
from damage rates measured at E r > E^. That is, to first order the average con- 
tribution of all thermal momentum directions to these damage rates is zero. 

Another effect that temperature may have on threshold energy concerns the problem 
of damage retention. Whenever thermal vibrations exist, there is a finite probability 
that the irradiation -induced interstitials will gain sufficient thermal energy to migrate 
from one lattice position to another until a vacancy is encountered and both defects are 
annihilated. Because of strain considerations, the activation energy and thus the tem- 
perature required for this migration should increase with interstitial -vacancy separa- 
tion, being smallest for those CPIA nearest their vacancies and largest for the FIA. In- 
deed, extensive low temperature irradiation studies on copper and other metals (ref. 17) 
have found that this model adequately explains the variety and behavior of the damage re- 
covery peaks observed during annealing through Stage I, the temperature region contain- 
ing initial damage recovery. That is, as temperature is increased within Stage I, CPIA 
with larger activation energies for migration gain sufficient thermal energy to perform 
correlated motion toward their own nearby vacancies (first-order kinetics); until at the 
end of Stage I the uncorrelated FIA wander freely throughout the lattice, encountering 
vacancies in a random fashion (second-order kinetics). Therefore, although both CPIA 
and FIA can be created at any irradiation temperature T, only those interstitials whose 
characteristic recovery temperatures are above T will be retained in the specimen to 
be measured. The only exception is the FIA which because of its random migration has 
a finite probability of being impurity-trapped and thereby retained above its intrinsic re- 
covery temperature T F j^. Generally Stage II, the next higher temperature region above 
TpiA» h as b een found to encompass those recovery peaks associated with FIA release 
from impurity traps. Therefore, if greater onset thresholds for larger interstitial- 
vacancy separations are assumed and TEDE is neglected, it follows that E^(T) should 
increase with increasing irradiation temperature within Stage I, reaching at T FF ^ the 
value E^(FIA) which should remain constant throughout Stage II. However, E^ as 
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determined with the full Stage II theory should be unaffected by the damage retention 
problem because, as explained with equation (4), the retention aspect is explicity con- 
tained in the probabilities f(e) and Jl - (2a^$/N u )J - */^ and not in the displacement 
probability cr d (e,E d ). 

In summary, the possible theoretical effects of a 350 K Stage II irradiation on thresh- 
old energy in tungsten are as follows. Onset threshold may or may not be larger than 
E°(0 K)-, that is, 

E°(350 K) = E°(0 K) + AE FIA - (2. 5±2. 5) = E d (FIA) - (2. 5±2. 5) (10) 

Here AE FIA is the difference in recoil energy required to create an FIA and that re- 
quired to create a Frenkel pair with minimum defect separation. The last term in paren- 
theses is the average contribution (eV) at 350 K of the TEDE. On the other hand, the ef- 
fective threshold E d as determined with the full theory (eq. (4)) should have the same 
value at any irradiation temperature within Stage II as that determined near 0 K (eq. (3)). 

Defect recovery in tungsten . - From the previous discussion it is obvious that 
proper analysis of the temperature dependence of threshold in tungsten depends on a 
knowledge of the temperature regions during which CPIA and FIA migrate and impurity- 
trapped interstitials are released. The initial resistivity recovery measurements on 
fast -neutron -irradiated tungsten by Kinchin and Thompson (ref. 18) have shown the exist- 
ence of three low temperature recovery stages: Stage I, 0 to 100 K; Stage II, 100 to 
600 K; and Stage in, 600 to 750 K. Thompson (ref. 1) suggested that Stage I recovery be 
assigned to interstitial motion, Stage II to the release of impurity -trapped interstitials, 
and Stage III to vacancy motion. Since then, further radiation damage experiments on 
tungsten have found general agreement with these assignments. 

Within Stage I, for example, the many recovery peaks observed by Coltman, 
Klabunde, and Redman (ref. 19) after thermal neutron irradiation are suggestive of the 
recombination of a variety of close Frenkel pairs each with different interstitial -vacancy 
separation. Regarding the location of T F j A , DiCarlo, Snead, and Goland (ref. 8) found 
that in electron-irradiated tungsten more damage recovered between 45 and 100 K after 
the specimen was predoped with excess vacancies. Such behavior would be expected if 
FIA move in this temperature range because the larger vacancy concentration should in- 
crease the probability that the FIA will be vacancy annihilated rather than impurity 
trapped. Neely, Keefer, and Sosin (ref. 20) have shown that recovery near 70 K in 
electron-irradiated tungsten followed second -order kinetics, a result which suggests FIA 
migration to vacancies near this temperature. More recently Scanlan, Styris, and 
Seidman (refs. 21 and 22) have performed a careful field ion microscope study with very 
high purity tungsten irradiated from 8 to 15 K with 20-keV tungsten ions. In the presence 
of the best imaging electric field they observed FIA reaching the specimen's surface 
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beginning at ~28 K and peaking at ~38 K. Because this surface field produced a large 
tensile stress within the bulk material, they were unable to conclude whether their HA 
annealing temperatures represented T FIA in a stress -free tungsten specimen. How- 
ever, they do point out that if the activation volume for FIA migration in tungsten were 
similar to that calculated by Johnson (ref. 23) for the ( 110) split FIA in a-iron (+0. 1 
atomic volume), then Tp^ in stress -free tungsten should occur near the 70 K recovery 
peak observed by Neely, Keefer, and Sosin. 

Support for the release of impurity -trapped interstitials in Stage II has been found by 
DiCar lo and Townsend (ref. 2) who observed internal friction peaks near 100 and 200 K 
in deuteron-irradiated tungsten. These peaks existed only in impure samples and were 
traceable to the stress -induced ordering of interstitial -impurity complexes. Regarding 
the mechanism responsible for Stage ni recovery, there exists some evidence (ref. 24) 
that a second form of the interstitial is annealing in this stage and that vacancy motion 
should be assigned to higher temperatures. Because the present results do not appar- 
ently give preference to either of these proposed mechanisms, we will not pursue this 
Stage III question further but only emphasize the generally accepted fact that tungsten 
vacancies do not migrate until at least 600 K or well above the maximum temperatures 
achieved here. 


Experimental Effects of Temperature 

Although there have been a number of low temperature radiation damage studies in 
tungsten, only two have yielded threshold values with accuracies comparable to that 
obtained here. Roberts (ref. 25) using electron irradiation at ~35 K measured 
damage rates near E^. Although he did not reach energies where no damage was 
observed, he was able to extrapolate his data to the zero damage point, obtaining 
E^(35 K) = 47±3 eV. Also, Neely, Keefer, and Sosin (ref. 20) measured electron- 
induced damage rates at three energies near 4 K, from which they estimate 
E°(4 K) = 40±2 eV. 

No previous experiment has determined a low temperature E^ for tungsten. Never- 
theless there seems to be sufficient published damage rate as a function of energy data to 
arrive at an approximate value for the effective threshold. For example, one can obtain 
the four data points shown in figure 3 from the 4 K damage rates of Neely, Keefer, and 
Sosin (ref. 20) and the 20 K rate of DiCar lo, Snead, and Goland (ref. 8) after correction 
for energy loss (eq. (6)) and path length (eqs. (7) and (8)). It should be emphasized that 
there can be significant deviations in electron dose measurements from laboratory to 
laboratory. In addition the 20 K point was probably observed above the annealing tem- 
peratures of some CPIA implying that, if the irradiation were performed near 4 K, a 
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larger damage rate would have been measured. To account for this latter problem we 
show by an arrow in figure 3 that this point could be higher near 4 K. Assuming a few 
values for this point and neglecting dose errors, we have fitted the three high energy 
damage rates to the theoretical displacement cross sections of Burke, Grosshard, and 
Lowe (ref. 14). The lowest energy damage rate was omitted because of the possibility 
of subthreshold displacements. We thus obtain the curves in figure 3, indicating 
E^ > 46 eV and p-y > 4. 1X10"° ohm -centimeter per atomic percent of Frenkel pairs. 

This latter result can be compared with the value determined by Schultz (ref. 7) of 
2. 5X10"° ohm -centimeter per atomic percent of tungsten vacancies. 

A summary of all results concerning the temperature dependence of threshold en- 
ergy in tungsten is given in table II. The defect recovery model which puts T FIA near 
70 K has been employed to label E^(35 K) as a Stage I measurement. As previously ex- 
plained with equation (5), the Stage II E^(350 K) may be in error due to the fact that it 
was determined under the assumption that f(e) was independent of energy. However, 
since E^ should be unaffected by temperature, the correct 350 K value which would 
have resulted from employing the full theory (eq. (4)) should be that determined near 4K. 
Thus it can be seen from table II that the assumption of constant f produced only a small 
error, if any, on the determination of E^ for tungsten. 

The majority of the data in table II indicates that across Stage I and a large portion 
of Stage II threshold energy in tungsten is nearly independent of temperature. The only 
exception appears to be the E^(4 K) result which is derived from the lowest energy point 
of figure 3. The fact that this point does not fall on the best cross section lines can be 
explained either by a large anistropy in E^hkZ) or by the existence of subthreshold dis- 
placements. For instance, if the threshold in one crystal direction is lower than the av- 
erage threshold for all directions, then near threshold a small amount of damage due to 
displacements in that specific direction will be observed; whereas at the higher energies 
the average threshold is exceeded and the damage is caused by displacements in essen- 
tially all directions. Or, the presence of "soft spots, " such as impurities and disloca- 
tions, may lower the binding forces surrounding the recoil atom, thus allowing displace- 
ments at smaller or subthreshold energies. This effect may also be due to displacement 
of impurities themselves (ref. 3). Considering the apparently strong effects of impuri- 
ties on defect recovery in tungsten (refs. 21 and 22) and the fact that tungsten's threshold 
is probably greater than that of any impurity it contains , one would tend to attribute the 
4 K lowest data point to impurity -related subthreshold displacements. Nonetheless, 
whatever its source, this type of low energy damage recovers somewhere between 4 and 
350 K since it was not observed in the results of figure 2. 
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Fundamental Aspects 


Having established a bare outline of the temperature dependence of threshold in tung- 
sten, let us consider what fundamental knowledge might be deduced by comparing the data 
of table II with displacement and recovery theory (see Theoretical Effects of Temperature 
section). 

Because the damage retained during a Stage n irradiation of tungsten and other 
metals can be attributed to equal concentrations of impurity -trapped FIA and vacancies, 
it follows that E^(350 K) = 50±2 eV is the minimum recoil energy required for FIA pro- 
duction at 350 K. Although the exact contribution of TEDE to this onset is unknown, the 
value of 50±2 eV should be the onset threshold measured at any Stage II temperature 
(100 to 600 K). This conclusion is based on the fact that damage retention anywhere 
within Stage II requires FIA production and that the only perturbation to onset is TEDE 
which changes insignificantly across Stage II (from -3 to -6 eV). 

To estimate interstitial thresholds in a tungsten lattice free of thermal vibration, 
one can employ equation (10) to find 

E d (FIA) = 53±5 eV (11) 

If it is assumed that the 4 K onset measured by Neely, Keefer, and Sosin was not caused 
by extrinsic sources, then 40±2 eV equals E^(0 K), the threshold for the production of a 
close Frenkel pair with minimum interstitial -vacancy separation. Thus, 

AE FIA = 13±7 eV (12) 

is the extra recoil energy required above E^(0 K) to deposit a tungsten interstitial in 
a lattice position at which its motion is unaffected by its own vacancy. This energy 
difference would be even smaller if E^(4 K) were due to subthreshold displacements. 

For if this were the case, one might then assume E^(0 K) = E^(35 K) (neglecting 
TEDE and the retention problem) to find AE^^ = 6±8 eV. 

Although only polycrystalline specimens were employed in the present study, general 
information concerning E^hkZ) for tungsten may be obtained by comparing the above re- 
sults with the computer calculations on bcc iron (ref. 16). For instance, let us assume 
that the following directional results of this latter study are also applicable to bcc 
tungsten: 

(1) E^(0 K) = E d ( 100) < E d (hkZ) ; that is, near 0 K the first observable damage is 
produced by atoms which have recoiled in the ( 100) directions. 

(2) E d ( 111) > E d (hkZ) ; that is, neglecting a few low probability recoil directions, 
atoms recoiling in the (111) directions require the largest recoil energy to produce an 
interstitial. 
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(3) E d (FIA) > E d (lll) ; that is, replacement collision chains in the close-packed 
(111) directions are the most probable source of FIA. 

Applying the previously mentioned iron results to the data of table n, one obtains 


and 


E d ( 100) = 40±2 eV < E d (hkZ) > 

► 

E d (hkZ) < E d (lll) < 53±5 eV^ 


( 13 ) 


indicating threshold energy in tungsten does not vary significantly with recoil direction. 
This variation would even be less if E d (4 K) were caused by subthreshold displacements. 
Approximate threshold isotropy is also suggested by the fact that one effective threshold 
fits the commercial data of figure 2 at all energies . If FIA threshold varied appreciably 
with recoil direction, E d over any small energy interval would not be independent of 
electron energy but would increase monotonieally as more directions for FIA displace- 
ment became active. 

The indication of approximate threshold isotropy in tungsten is in opposition to the 
theoretical estimates of Andersen and Sigmund (ref. 26) who calculate 

E d ( 111) =2.4E d (100). Admittedly the directional dependence of E d is best determined 
with single crystals where directional measurements are more straightforward and are 
not subjected to the uncertainties contained in the association of certain defect types with 
certain crystal directions. Nonetheless, the previous analysis does suggest that 
E d (hkZ) in tungsten is more isotropic than the calculations of Andersen and Sigmund 
would seem to indicate. 


DAMAGE PRODUCTION ABOVE THRESHOLD 

According to low dose Stage n theory (eq. (5)), the energy dependence of the 350 K 
damage rate is produced by two factors: f(e) and cr d (e,E d ). Since E d has been shown 
to be temperature independent, it follows from equation (3) that cr d has the same energy 
dependence as the damage rate near 4 K. Thus f(e), the FIA fraction of the total inter- 
stitial concentration, is directly proportional to the ratio of the 350 K to the 4 K rates. 
Explicity from equations (3) and (5), 
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ajc 

Applying this equation with p = pp to the data of figure 3 and the commercial results 
of figure 2, one can then determine the energy dependence of f(e). For instance, at 
1. 6 MeV, f = 0; at 1. 9 MeV, f = 0. 22; and at 2.4 MeV, f is 0. 33 or less depending on 
whether is 46 eV or greater. Thus, as energy increases above threshold, f(e) 
rises rapidly from zero and then tends to level off. Since f is the fraction of FIA, it 
follows that after a low temperature irradiation the percentage of Stage I recovery due to 
FIA should also show a tendency to saturate with increasing energy. Energy dependence 
such as this has already been observed for FIA recovery in some fee metals (ref. 27). 

The 350 K damage rate results indicate that f may also depend on impurity concen- 
tration. Assuming that the commercial (C) and high purity (HP) specimens have the 
same and p* , then one finds from low dose theory (eq. (5)) and figure 2 that 
f Hp/ f C — 0.5. Because neither specimen type should have contained enough impurities 
to affect the displacement process above threshold (ref. 3), this purity effect must be 
related to the process of FIA retention. Accepting this, the simplest interpretation 
which can then be invoked is that there is a higher probability with larger impurity con- 
centration that some CPIA will be created within the capture radii of the impurities. Up- 
on creation in Stage II these CPIA are immediately attracted away from their own vacan- 
cies by the impurities and thereby counted as trapped FIA. This mechanism is consist- 
ent with observed impurity effects in copper (ref. 28) where close-pair recovery in 
Stage I has been measurably reduced by large impurity concentrations. However, even 
though this CPIA capture mechanism is a reasonable interpretation, there are at least 
two other possible sources for the purity effects. These are based on the possibility that 
the inequality in damage rates for the two specimen types was produced totally, or in 
part, by the parameter p* . For instance, due to a difference in impurity types the re- 
sistivity of a trapped interstitial in the pure specimen may be smaller than that in the 
commercial specimens. Or, because defect resistivity was measured at 291 K, devia- 
tions from Mathiessen's rule may exist, causing > p^p. Nevertheless, whatever 
the source, it is apparent that in the low dose region resistance ratios ranging from 10 
to 1000 produce only a small effect on Stage n damage rates. 

Finally, to determine the energy dependence of low dose Stage H defect production 
in tungsten, one can employ the following relation which adequately describes the results 
of figure 2: 



( 15 ) 


^ (e) =f oP *a d (e,52eV) 

Ld$ J n 


where the product f Q p* is a constant. If the energy dependence of f is averaged from 
1. 9 to 2.4 MeV, then equation (15) gives f^ = 0. 28±0. 06 and fjjp = 0. 14±0. 03. Here it 
is assumed that 46 eV < E*[ < 52 eV, p* = p F , and Mathiessen’s rule holds at 291 K. 
Thus in the low dose region 


L^(e) =J_ dp(e)l = (0.22±0.12We,52 eV) (16) 

L d# J n p* Ld$ J n 

where N F is the electron-induced Frenkel pair concentration created near 350 K in 
polycrystalline tungsten of any practical purity. If variations in p* are neglected, the 
contents of the section Theoretical Effects of Temperature predict that equation (16) 
should hold across all of Stage II, that is, 100 to 600 K. Further experimentation within 
this temperature range should test this prediction. 


THRESHOLD IN bcc TRANSITION METALS 

The near isotropy of E d (hkZ) in elastically isotropic tungsten suggests that perhaps 
some correlation exists between threshold energy and directionally dependent elastic 
properties. Although there seems to be no a priori reason why such a relation should 
exist, it was of interest to check whether similar behavior could be found in other bcc 
transition metals. Lucasson and Walker (ref. 29) who electron-irradiated both polycrys- 
talline iron and molybdenum specimens near 20 K found that E^ for iron was not well- 
defined; that is, the energy-dependent data could not be fit to a simple step-function prob- 
ability of displacement but required a function which rose from zero to unity over a wide 
range of recoil energies. Erginsoy, Vineyard, and Englert (ref. 16) subsequently 
showed that this requirement of a nonsimple probability function could be explained by 
anistropy in iron’s E d (hkZ). Regarding molybdenum, however, Lucasson and Walker 
found that a simple step-function displacement probability could fit all the data, even 
those damage rate points close to zero. Thus by comparison with the iron it appears 
that threshold for molybdenum is isotropic; that is, E d (hkZ) = E^. Because iron and 
molybdenum are elastically anisotropic and isotropic, respectively, the results of 
Lucasson and Walker seem to qualitatively support a correlation between the directional 
dependence of threshold energy and elasticity. Encouraged by this we have examined 
whether any general quantitative agreement could be found among the three bcc metals 
for which thresholds have been measured. 


17 



An elastic property which like threshold depends on crystal direction is Young’s 
modulus Y. For the stress direction (hkZ) in a cubic crystal with elastic constants 

C ll’ C 12’ and C 44’ Y is ( ref - 3 °) 

(Y(hkZ) ) -1 = (Y(IOO))" 1 - (A - l)r(C 44 ) -1 (17) 


Here A, the anisotropy factor, equals 2C 44 /(C^ 1 - C 12 ); 


Y(100) = 


( c n - c i 2 )( c ii + 2 C i 2 ) 

( c n + C 12 ) 


(18) 


and r, the orientation factor, is expressed in terms of the direction cosines y^, y 2 , 
and y 2 of the stress direction relative to the three cube axes 

_ 2 2 2 2 2 2 

r = YjY 2 + r 2 r 3 + y 3 y i (19) 


An anisotropy factor of unity indicates elastic isotropy, that is, Y(hkZ) = Y(100) . Using 
the appropriate low temperature elastic constants (ref. 31), we calculated Y(100)(F = 0) 
and Y( 111) (r = 1/3) for a number of bcc transition metals. The results are given in 
table in Also included are the available experimental for each metal plus the 
computer -calculated E d (hkZ) for iron (ref 16). The uncertainty in E d ( 100) for tung- 
sten arises from the possibility that subthreshold displacements may have influenced this 
measurement. Comparing the experimental E^ with the minimum Y(hkZ) we find 

jjjO 

— — s 1. 2X10" 6 (eV)(cm 2 )/N = 20 A 3 (20) 

Y . 
min 


Assuming this threshold to modulus ratio holds for all bcc transition metals and all crys- 
tal directions, we can then predict E^hkZ) from Y(hkZ). These empirical predictions 
are given in the second last column of table III. Finally, the last column indicates the 
Andersen and Sigmund (ref. 26) theoretical estimates of E^hkZ) for each metal. 

From table HI it can be seen that wherever experimental or computer -calculated bcc 
thresholds are available, fairly good quantitative agreement is achieved by assuming 
E^hkZ) directly proportional to Y(hkZ) through the equation (20) constant. It is also 
clear that this type of empirical correlation is more useful than the predictions employed 
by Andersen and Sigmund. These authors assumed that each bcc E^(hkZ) is propor- 
tional to the value of the interatomic potential at an atomic separation equal to the near- 
est neighbor distance along the (hkZ) direction. Their directionally dependent proportion- 
ality constants were obtained from the computer results on iron; whereas their 
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interatomic potentials were assumed to be Born-Mayer, that is, A Q exp(-r/a Q ). Regard- 
ing the parameters a Q and A Q , they found that good agreement with the threshold re- 
sults of fee copper, silver, and gold could be achieved by assuming a Q is a constant 
and A q is proportional to the 3/2 power of the atomic number. From table HI it is ap- 
parent that this choice of parameters for the bcc transition metals is not as universal as 
Andersen and Sigmund assumed. 

Although reasons for the apparent success of threshold -modulus correlation are far 
from obvious, perhaps the following argument may shed some light. In general, appli- 
cation of correlation with one proportionality constant to the fee metals has proven to be 
less successful than the method of Andersen and Sigmund. One possible exception is 
nickel which, although it possesses an fee structure, is a transition metal. This sug- 
gests that the unfilled d-shells are the source of threshold -modulus correlation. In any 
metal, threshold should be the energy required to force the outer shells of the struck 
atom through the outer shells of its nearest neighbors. In the transition metals, cohe- 
sive energy and thus Young's modulus are due primarily to the local interaction of the 
unfilled d-shells. However, in the nontransition metals such as the noble metals, cohe- 
sive energy and modulus are due to the nonlocal long range interactions of the valence 
electrons. Therefore, threshold -modulus correlation might be expected in the transi- 
tion metals but not in the noble metals. This possibly might also explain why the Born- 
Mayer constants derived by Andersen and Sigmund from the noble metals do not apply 
as well to the bcc transition metals. Suffice, to say, the correlation method like the 
Andersen -Sigmund method is based on only three metals so that future bcc threshold 
measurements may totally nullify its general usefulness. Regarding this point it appears 
that a good test metal would be chromium since the E^j prediction by modulus correla- 
tion (31 eV) is noticeably different than that given by Andersen and Sigmund (12 eV). 


CONCLUSIONS 


The present study has measured the energy dependence of the rate at which electron- 
induced radiation damage accumulates in commercial and high purity tungsten near 350 K. 
From these results the minimum or threshold energy for the onset of observable damage 
was determined to be 1. 62±0. 04 MeV, or, in terms of atomic recoil energy, E~j = 

50±2 eV. An "effective" threshold of 52±2 eV was also determined by directly compar- 
ing the energy dependence of the damage rates to the energy dependence of the "best fit" 
theoretical cross section calculated from a step-function displacement probability. 

The aforementioned results in combination with low temperature damage rate data, 
bcc displacement theory, and the accepted mechanisms for defect recovery leads one to 
the following conclusions: 


19 



1. Since the irradiations were performed at temperatures within Stage II recovery, 
the observed damage consisted primarily of equal concentrations of vacancies and 
impurity -trapped interstitials. In order to reach impurities, these interstitials were 
initially created at lattice positions from which they could migrate freely, that is, unaf- 
fected by their own vacancies. Thus the onset threshold result of 50±2 eV is attributed 
to the minimum recoil energy required to produce an uncorrelated or "free" interstitial 
atom (FIA) in tungsten near 350 K. Because the effect of thermal vibrations varies in- 
significantly across Stage II, 50±2 eV should be the onset threshold measured anywhere 
from 100 to 600 K. 

2. From a re-examination of previously published low temperature damage rates, 
the effective threshold near 4 K was determined to be 46 eV or greater. This result 
suggests that the 4 K onset of 40±2 eV measured by Neely, Keefer, and Sosin may have 
been caused by impurity -related subthreshold displacements. 

3. For a tungsten lattice free of thermal vibration, the minimum recoil energy to 
create an FIA is 53±5 eV. Comparison of this result with E^j(4 K) indicates that the 
extra recoil energy AEpj^ = 13±7 eV is required to increase the minimum interstitial- 
vacancy separation to the point where the interstitial is uncorrelated or unaffected by the 
vacancy. If subthreshold events did affect the measured 4 K onset, this energy differ- 
ence is even smaller, that is, AEp-^ = 6±8 eV. By assuming that the qualitative com- 
puter findings on bcc iron apply also to bcc tungsten, it can be shown that the variation 
in E d (hkZ) with recoil direction is no greater than AE FIA . Thus threshold in tungsten 
is nearly isotropic. 

4. A difference of two orders of magnitude in impurity content was found to have 
only a small effect on low dose damage rates. Consequently, an equation is derived for 
estimating the Stage II energy dependence of Frenkel pair production in tungsten of any 
practical purity. This result, however, is limited to low doses and recoil energies 
below ~100 eV. 

5. Finally, the finding of approximate threshold isotropy in elastically isotropic 
tungsten has suggested a possible correlation in bcc transition metals between the recoil 
direction dependence of threshold energy (E^hkZ)) and the directional dependence of 
Young's modulus (Y(hkZ) ). Indeed, it is shown that both qualitative and quantitative 
evidence exist for such correlation in iron and molybdenum, the only other bcc metals 
for which threshold has been determined. The low temperature directionally dependent 
threshold energies of iron, molybdenum, tungsten, niobium, chromium, tantalum, and 
vanadium are calculated using the relation E^(hkZ) /Y(hkZ) = 1. 2x10”® (eV)(cm^)/N. 
These results are compared with the theoretical predictions of Andersen and Sigmund. 
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Wherever experimental data exist, significantly better agreement is found for the corre- 
lation method. Only further experimental and theoretical studies of bcc thresholds will 
determine the general usefulness of the modulus -correlation approach. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 3, 1971, 

129-03. 
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TABLE I. - ENERGY DEPENDENCE OE DAMAGE RATE NEAR 350 K 


Specimen 

Incident electron 
energy, 

e i- 

MeV 

Maximum recoil 
energy , 

E max 
r ’ 

eV 

Damage rate, 
(dp/d<l>) m , 

ohm-cm/electrons/cm* 2 

Commercial 
purity 1 

2. 30±0. 01 

89 

a (7.83±0.39)xl0' 27 

Commercial 

1. 70±0. 01 

54 

(0. 34±0.04)xl0' 27 

purity 2 

1. 90 

65 

(1. 40±0. 14) 


2. 10 

77 

(4 . 00±0. 15) 


2. 20 

83 

(5. 22±0. 14) 

Commercial 

1. 59±0. 01 

49 

( 0 . 00 ± 0 . 01 )xl 0 ' 27 

purity 3 

1.85 

62 

(0. 92±0. 03) 


2.00 

71 

(2. 69±0. 04) 


2. 25 

86 

(5. 68±0. 11) 


2.40 

96 

(8. 11±0.01) 

High purity 1 

2. 10±0. 01 

77 

(1.59±0.03)xl0~ 27 


2.40 

96 

(4. 52±0. 24) 


a No temperature controller. 


TABLE II. - TEMPERATURE DEPENDENCE OF 
THRESHOLD ENERGY IN TUNGSTEN 


Threshold energy 

Stage I 
(~4K) 

Stage I 
(—35 K) 

Stage II 
(-350 K) 

Threshold for onset of damage, 

a 40±2 

b 47±3 

c 50±2 

E° eV 




d’ 

Effective threshold, 

d >46 


e 52±2 

E d' eV 





a Ref. 20. 
b Ref. 25. 
c Present study. 

d From combination of refs. 20 and 8; see text and fig. 3. 
e Present study under assumption f(e) = constant. 





TABLE m. - COMPARISON OF THRESHOLD ENERGIES IN bcc TRANSITION METALS 


Metal 

Crystal direction, 
(hkZ) 

Young's modulus 
(experimental) a , 
Y(hkZ), 
N/cm 2 

Threshold energy, E^hkZ) , eV 

Experimental 

(approximate) 

Theoretical 

Present 

study 

Andersen- 

Sigmund® 

Iron 

(100) 

14x10® 

c 17 

17 

14 


(110) 

23 

d 34 

28 

-- 


(111) 

29 

d 38 

35 

31 

Molybdenum 

(111) 

32X10 6 

e 37 

38 

35 


(100) 

35 

e 37 

42 

15 

Tungsten 

(100) 

42x10® 

f 40 to 47 

50 

34 


(111) 

42 

g 53 

50 

80 

Niobium 

(111) 

8xl0 6 

— 

10 

26 


(100) 

15 

— 

18 

11 

Chromium 

(111) 

26x10® 

— 

31 

26 


(100) 

36 

— 

43 

12 

Tantalum 

(100) 

15X10° 

— 

18 

25 


(111) 

23 

... 

28 

61 

Vanadium 

(111) 

13X10® 

— 

15 

-- 


(100) 

15 

— 

18 

— 


a From eqs. (17) -(19), and ref. 31. 
b Ref. 26. 
c Ref. 3. 

^Computer -calculated results from ref. 16. 
e Deduced from ref. 29; see text. 
f Refs . 20 and 25. 

^Present study. 




Incident 

electron 



Figure 1. - Typical damage production curves for low dose 
irradiations near 350 K. Commercial purity specimen 3. 



Figure 2. - Damage rate near 350 K plotted against average 
electron energy for 2-mil commercial and 5-mil high 
purity tungsten wires. All data have been corrected for 
path length (eqs. (7) and (8)). 



Figure 3. - Damage rate near 4 K plotted against average electron 
energy as measured by Neely, Keefer, and Sosin (ref. 20) with 
3-mil tungsten wires. Damage rate near 20 K measured by 
DiCarlo, Snead, and Goland (ref. 8) with 5-mil wire is also in- 
cluded for comparison. All data have been corrected for path 
length (eqs. (7) and (8)). 
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